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Abstract 

We consider the minor process of (Hermitian) matrix diffusions 
with constant diagonal drifts. At any given time, this process is de- 
terminantal and we provide an explicit expression for its correlation 
kernel. This is a measure on the Gelfand-Tsetlin pattern that also 
appears in a generalization of Warren's process [26], in which Brow- 
nian motions have level-dependent drifts. Finally, we show that this 
process arises in a diffusion scaling limit from an interacting particle 
system in the anisotropic KPZ class in 2 + 1 dimensions introduced 

mm. 



1 Introduction 

In this paper we determine a determinantal point process living on the 
Gelfand-Tsethn cone GT^v, 

GT^ = {{x\ x^,...,x^) eR' xR^ X ■■■ xR^ : x^^' < < xlX\} (1) 

arising both from random matrices diffusions and and interacting particle 
systems. An element x G GTat is called a Gelfand-Tsetlin pattern. Here 
is a graphical representation of a; G GT4, which illustrates the interlacing 
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Measures on Gelfand-Tsetlin patterns naturally appear in several fields of 
mathematics like (a) random matrix theory [2], [181 IHl EH] where the ques- 
tion of universality was recently approached in [20], (b) random tiling prob- 
lems [3 [221 [22], (c) representation theory [H1[H], and (d) interacting particle 
systems [H, [211 [22] and diffusions [2S1 [21] • Probably the most famous exam- 
ple which belongs to more than one of these classes is the Aztec diamond. 
Indeed, the measure on the Aztec diamond both comes from a random tiling 
problem [11] and can be obtained through a Markov chain [21] on a (discrete) 
Gelfand-Tsetlin pattern, which itself can be seen as a special case of the more 
general Markov chain construction in |1]. A continuous space analogue is for 
instance Warren's process, a system of interacting standard Brownian dif- 
fusions [261 EZ]- Recently, the link between the "drift-less" case of [1] and 
Warren's process has been studied in [T5] . 

In this paper we consider CUE matrix diffusions with drifts. Its eigen- 
value process at a fixed time is a determinantal point process and we explic- 
itly determine its correlation kernel. We then show that the point process 
arises from a diffusion scaling limit of an interacting particle system in the 
anisotropic KPZ class in 2-1-1 dimensions [5l[l], as well as in a generalization of 
Warren's process [2S] if we let the Brownian motions to have level-dependents 
drifts. 

For the zero-drift case, the GUE minor process was described in [18], 
while Warren's process was introduced in [26], and it is known that the two 
processes coincide along "space-like" paths [12], they both are Markovian 
and determinantal, but in the whole "space-time" the two processes are dif- 
ferent [T]. Here we focus on the fixed-time process although the connection 
will certainly hold along "space-like" paths as for the zero-drift case. 



Matrix diffusions 

The first model we study on GTn is a variant of GUE minor process which 
has been introduced in [18]. Consider an x iV Hermitian matrix H with 
eigenvalues < ■ ■ ■ < Aj^. Denote by the submatrix obtained by 
keeping the first n rows and columns of H, and its ordered eigenvalues by 
A" < ■■■ < A^. The collection of all these eigenvalues (A^,...,A^) then 
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forms a Gelfand-Tsetlin pattern, with A" = (A", . . . , A"). In this paper we 
take H{t) to be a GUE matrix diffusion perturbed by a deterministic drift 
matrix M = diag(/ii, . . . ,fiN), i-e., we consider G{t) = H(t) + tM with H 
evolving as standard GUE Dyson's Brownian Motion starting from 0. The 
eigenvalues point process ^ has support on M x {!,..., A^} and it is given by 

^(dx,m)= ^ Sn,mSx]^{dx) (2) 

l<k<n<N 

and their correlation functions are given as follows (for the proof, see Sec- 
tion [2]). 

Theorem 1. For a fixed time t > consider the eigenvalues' point process 
on the N suhmatrices of H[t). Then, its m-point correlation function is 
given by 

^"((xi, ni), . . . , {xm, rim)) = det[Kt{{xi, Ui), {xj, nj))]i<ij<m, (3) 
with {xi, Uj) G M X {!,..., A^} and correlation kernel 

n' 

UM, {x',n')) = + K-,{x)<l>iUx'), (4) 

k=l 

wher^ 

(T\n'~n p z{x'-x) 

(^("'"')(x,xO= ^ ^ . / dz- ^ ^ ![„<„,], (5) 



wzi/i /i_ < min{/ii, . . . , /iat}. 

Remark 1.1. The integral for (j)'^'^'^'^ in ([5]) is only well-defined forn' —n > 1. 
For n' — n = 1 we set x') := = e^"*^^ ~^^l[x>x'] instead. 



^ For a set S, the notation ^p_^ dw f{w) means that the integral is taken over any 
positively oriented simple contour that encloses only the poles of / belonging to S. 
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Warren's process with drifts 

Our second model is Warren's process with drift which describes the dynam- 
ics of a system of Brownian motions {B^, 1 < k < n < N} on GTtv, where 
Bl is a standard Brownian motion with drift fii starting from the origin. The 
Brownian motions Bf and 5| are Brownian motions with drift H2 conditioned 
to start at the origin and, whenever they touch Bj, they are reflected off Bl. 
Similarly for n > 2, BJ^ is a Brownian motion with drift /i„ conditioned to 
start at the origin and being reflected off by -B^~^ (for k < n — 1) and B^zl 
(for k > 2). The process with /ii = ■ ■ ■ = /xat = was introduced and studied 
by Warren in [26]. In this way we obtain a system of Brownian motions with 
drifts that lies in GTjv- The correlation functions of this process at a flxed 
time agree with those of the perturbed GUE minor process (the proof is in 
Section [3]). 

Theorem 2. For a fixed time t > consider the point process of the positions 
of the Brownian motions {BJ^it) : 1 < k < n < N} described above. Then, 
its m-point correlation function is also given by ^ . 

Interacting particle system 

Finally we introduce a discrete model which gives rise to Warren's process 
with drifts under a diffusion scaling limit. This model is a generalization 
of TASEP with particle-dependent jump rates [S] to the 2-1-1 dimensional 
particle system with Markov dynamics introduced in [1]. We denote by 

e Zi the position of a particle labeled by (/c, n), with 1 < k < n < N , and 
call n the "level" of the particle. Particle (/c, n) performs a continuous time 
random walk with one-sided jumps (to the right) and with rate v^- Particles 
with smaller level evolve independently from the ones with higher level like in 
the Brownian motion model described above. More precisely, the interaction 
between levels is the following: (a) if particle (/c,n) tries to jump to x and 
x'^Zi = X, then the jump is suppressed, and (b) when particle {k,n) jumps 
from X — 1 to X, then all particles labeled by {k + i,n + t) (for some £ > 1) 
which were at x — 1 are forced to jump to x too. This is a particle system 
with state space in a discrete Gelfand-Tsetlin pattern. 

Consider the diffusion scaling with appropriate scaled jump rates 

t = rT, xl = rT-^\l, ^„ = 1 _ (8) 

Then, in the T ^ 00 limit, the particle process {x'^{t)} converges to the 
GUE minor process with drift {A^(r)}. 
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More precisely, let us denote by P^' the probability measure on these 
particles with jump rates v = {vi, . . . , vjy) given in ([H]). We fix r > and set 

z/t(A) = r /_ ^fc(^^-^^ G for all 1 < A; < n < Nj (9) 



where C M are Borel sets, A = ni<fc<n<Af^fc- Moreover, we define 

u{A) = P^(A^(r) G for all 1 < A; < n < A^) (10) 

where is the GUE minor measure with drift diag(/xi, . . . , /iat). In Section H] 
we show the following result. 

Theorem 3. As T ^ oo, uj- converges to v in total variation, i.e., 

lim sup \ut{A) - iy{A)\ = 0. (11) 



T-5.00 



JV(]V+l)/2 



A Borel 

In particular, — )■ z/ weakly. 
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2 GUE minor process with drift 
2.1 Model and measure 

Let {H{t) : t > 0) be a process on the N x N Hermitian matrices defined by 

(bkk{t) + fikt, iil<k<N, 
Hkiit) = I -^{bkeit) + ihe{t)), ifl<k<£<N, (12) 
[j=,{bUt)-ihi{t)), ifl <i<k<N, 

where {bkk, bke, bke} are independent one-dimensional standard Brownian mo- 
tion^. Denote by M = diag(/ii, . . . ,/iAr) the diagonal drift matrix added to 
the matrix H. Then, the probability measure on these matrices at time t is 
given by 

F{H G dH) = const x exp f-^^^^^i-J^^^ dH (13) 



^Herc, standard Brownian motions start from and arc normalized to have variance t 
at time t. 
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where dH = YYl^ dH,,Yli<, 

<k<N ^^^i'^j,k)dli^{Hjk) and const is the nor- 

mahzation constant. 

Since we are interested in the statistics of the eigenvalues' minors at time 
t, we first determine the measure on the eigenvalues of the N x N matrix. 

Lemma 2.1. Assume that /ii, . . . , fij\f are all distinct. Then under fll3p . the 
joint probability measure of the eigenvalues Ai, . . . , A^r of H is given by 

P(Ai e dAi,...A^ e dA^) 

= const x(f] e-(A-*M.)V(2t)^ ^i^^u ■ ■ ■ , >^n) ^ ^^^^ 

with const a normalization constant and A(xi, . . . ,Xm) = ni<i<j<m(^i 
the Vandermonde determinant. 

Remark 2.2. If fii, . . . ,fiN are not all distinct, we have to take limits in 
f|T4|) . For instance, if /ii = ■ ■ ■ = /iAr = /i, then 

P(Ai e dAi,...A^ e dA^) 

= const X (fl e-(^^-*'')'/(2*)') A2(Ai, . . . , A^) dAi ■ ■ ■ dA^. (15) 

Proof of Lemma \2.1[ We diagonalize H = U AU* with a unitary matrix U 
and the diagonal matrix A = diag(Ai, . . . , Aat). Then, 

^_Tr{H-tAf)V(2t)^^ = const X g" (A) df/ dA, (16) 

where dU is the Haar measure on the unitary group U. Moreover, since 

Ti{UAU* -tM) = TrA2 + t2TrM2-2tTr(f/Af/*M) (17) 

by integrating over U in f|T6|) and using the Harish-Chandra-Itzykson-Zuber 
formula we obtain the desired expression. □ 

Now we focus on the minor process. For 1 < n < let us denote by 
H^{t) the nx n principal submatrix of H{t) which is obtained from H{t) by 
keeping only the n first rows and columns. In particular, H^{t) = i?ii(t) and 
H^{t) = Hit). We denote by A^(t) < ■ ■ ■ < Xlit) the ordered eigenvalues 
of H^{t). It is then a classical fact of linear algebra that at any time t, the 
process (A^, . . . , A^)(t) lies in the Gelfand-Tsetlin cone of order A^, 

GT^ = {(a;\ . . . , x^) G X ■ ■ ■ X : x" ^ for all 1 < A; < A^ - 1}, 

(18) 
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where x"' ^ x'^'^^ means that and x'^'^^ interlace, i.e., 

<Xk< Xktl for all 1 < A; < n. (19) 
The induced measure on{A^: l<A;<r2< A^} is the following. 



Proposition 2.3. Fixt > 0. Then, under the measure f[T3l) . the joint density 
of the eigenvalues of {H^ : 1 < n < N} on GTjv is given by 



N N 



constx J]e-*'^^/2-Qg-(A-)V(2<)^(;,^) "Q ^m^a^ "Q ^-,^K-\ (20) 

fe=l k=l l<n<N 2<n<N 

l<k<n l<fc<n-l 

where the normalization constant does not depend on /ii, . . . ,fiN- 

Proof. We first derive (l20i) under the assumption that the /ii, at are all 
distinct. The case where some of the fii are equal is then recovered by taking 
the limit. We prove the statement inductively and follow the presentation in 
[13]. For = 1, the density is clearly proportial to exp(— (Aj — nitY /{2t)). 
For iV > 2, we consider an x matrix distributed according to f|T3|) 
which we write as 



where M'^~^ denotes the (A^ — 1) x (A^ — 1) principal submatrix of M, 
w G C^~^ is a Gaussian vector and x G M is a Gaussian variable. Then we di- 
agonalize H^~^, i.e., we choose a unitary matrix U such that H^~^ = UAU* 
with A = diag(A{^~^, . . . , A^lJ) the diagonal matrix for the eigenvalues. Since 
the Gaussian distribution is invariant under unitary rotations and w is inde- 
pendent of we have 



where = denotes equality in distribution. Applying the map ^ i— ?■ (A, f/), 
the get that measure f|T3|) on is proportional to 



exp I Tr 



W* X I ^ \ llN 



A2(A^-i) 
X dt/du;da;dA^-\ (23) 



7 



where dU is the Haar measure on the unitary group IAn^i. We consider 
only the part of fl2^ that depends on U and integrate over U^^i using the 
Harish-Chandra-Itzykson-Zuber formula 

f df^^TKA^.M--^/) ^ ^^^^^ ^ /^fyr^'^'^'- (24) 

A(A^-i)A(/.i,...,/.^_i) 
After this integration the measure (1231) reads 



N-l 



const X P(A^^i G dA^-i) e"^^-*'^^/^ "Q Q-\-'k\Vt^^ (25) 



fe=i 



We focus on the measure on Wk and represent the variables in polar coordi- 
nates, Wk = r^e^'^* with G M_|_ and ^ [0, Stt). Since the Jacobian of this 
transformation is given by ri ■ ■ ■ rjy-i, we get 



N-l N~l 



W e-\^>'\"l'dwk = n ^fce"^^/*dr, d^k: (26) 



fc=i fc=i 



where dr^ and dyp^ are Lebesgue measures on M+ and [0, 27r). Then we can 
express and a; in terms of the eigenvalues of H^~^ and if^, see e. g. 
for details, 

nf.,(Ar'-A?') 

if 



'FPV^I Tw^n ^lAA^-MA^^], 



TV 7V-1 

x = Tr(if^-i7^-) = 5:Af-5:Ar. 

4=1 k=l 

The Jacobian of the transformation T : (ri, . . . , r^r-i, x) t— )■ A^ is then given 
by 



A(A^^ 
A(A^ 



ri---r7v-i|detT'| = ^j—^I^^n-i^^n^, (28) 



and hence, given A^ ^, we have 



N-l N N-l 

fc=l fc=l fe=l 



X 



A(A 



N^ 



:^l[Aiv-M,iV]dA^M<^. (29) 
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Here we used that 2(r^ H h J = Tr(if^)2 - TT{H^-^f. Moreover, by 

the induction assumption for — 1 we have 



N~l N-l 



G dA^-1) = const X II e-'^"/' J] e-(^^ )'/(2t)^(;^Ar-i) 

fc=l k=l 

N-l 

X JJ e^"^^ Yl e-^"^r' JJ dA". (30) 



l<n<Af-l 2<n<Ar-l n=l 

l<A,<n l<fc<n 

Finally, inserting fl29p and (13 Op into (I25p and integrating out (f (which multi- 
plies the measure by a finite constant) results in the claimed result fl20|) . □ 

2.2 Correlation functions 

Now we want to determine the correlation functions of the point process 
on the eigenvalues {A^:l<A;<n< N}. We proceed by writ- 
ing the density in (^U^ as a product of determinants. Therefore we set 
(pnix.y) = e^"-^y~^^l{x>y} and introduce "virtual" variables A""^ = virt with 
the property that 0„(virt,?/) = e^"^. Then in ( 120|) we have, up to a set 
measure zero, 

n n—l 

det[MK~\>^l)]l<^,J<n = (31) 

j=i i=i 

Moreover, for A; = 1, . . . , we set 



where p„ are symmetric polynomials of degree n in n variables defined by 
Po = I and 

(—1)" f 2 

Pn{Xi, ■ ■ ■ , Xn) = . / dty e*" ''^(tf — Xi) • ■ ■ (tf — X„) for 72 > 1. (33) 



iv27r Jm 
Hence we have that 

N 

n e-(^^)V(-) A(A^) = const x det [K'-.m ,<.,<^, (34) 
fc=i 

which means that we can rewrite fl20|) as 

TV N 

const X n det [0„(Ar \ A,")] J] ^et i<,,<^- (35) 

n=l k=l 

9 



Note that by a change of variable w = {tz — x) / \/t we have 

(_T\N-k r 

K'U^) = / dze'^'/'-^^{z - ■■■{z- /x^). (36) 



A measure of the form fl35|) has determinantal correlation functions and the 
kernel can be computed with Lemma 3.4 of [6], see the appendix. 

2.3 Proof of Theorem [1] 

We prove the theorem first for /ii < ■ ■ ■ < /ijv and then use analytic contin- 
uation. Note that for n = N, the function ^I/"!!^ in ([6]) is the same as ^1/^1^ 
in ([36]). 

Lemma 2.4. The following identities hold. 

(i) For all n e {I, ... , N], k eZ andt > 0, we have * ^^i^^ = "^"^^Itk- 
(a) For n < n' , we have 4>n+i * ■ ■ ■ * 0n' = with given in 

Proof. Because of Re ^ < Hn we can exchange the two integrals, 

(0n*C:^,)(x) 

= r dye^-^y-^^^-^^ [ aze'^'/'-y^ ['~^'l---['~^''l 

J~oo 27ri [Z - fli) . . . [Z - fik) 

2vri iiM+^„ (2; - /ii) ... (2 - /ifc) 

-1)"^^"'= r , ,t.V2-x.. (2: - /il) ■ ■ ■ (^ - /in-l) 



I dze 

2vri iiM+^_ (2; - /xi) . . . (z - /ifc) 

This proves the first statement. To show (ii), we first consider the case 
n' — n = 2. A simple calculation gives 

( ll„{x'-x) lJ,n-l{x'-x)\ 

(38) 

which has the following contour integral representation, 

0("^2'")(a;,a;') = — / dz — — -. (39) 

27ri iiK+u_ [z - /i„_l) (z - /in) 
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For n' — n > 2, we get inductively that 
' ' ' d?/e^"(^-") / 



27ri i„oo iiR+u_ (2; - ^„+i) ■ ■ ■ (2; - 



d^- ^ / dye^^'^"-") (40) 



2vri [z - iin+i) ■ ■ ■ [z - llr. 

d2; 



27ri 7iR+;,_ (2; - - /i„+i) ■ ■ ■ (2; - /i„/) 

= (/.('^-i'"')(x,x'), 

where, as before, we could exchange the integrals because of Re 2; < /i„. □ 

Next we consider the n-dimensional space spanned by the set of func- 
tions 

{(/.I * (X?, ■),..., 0n-l * </.("-^'")(xr?, ■), ■)}• (41) 

According to Lemma 3.4 of [6] we need to find a basis {^"^1^ : 1 < A; < n} of 
that is biorthogonal to the set {^I/^il^ '-^ < n}, i.e.. 



/dxC:!fc(:^)C-£(^) = 5fc^, l<A:,£<n. 



(42) 



The form of the biorthogonal functions can be, with some experience, guessed 
from the form of the kernel [TU] . 

Lemma 2.5. We have: 

(i) Vn is spanned by {x h-> e'^*^ : 1 < k <n}. 

(a) The functions {$"1!^ ■.l<k<n} are given by ([7]). 

Proof. For any e > we have 

(0,*0('=''^))(x^\x) 



dye'^^y^ / dz- ^ -. (43) 



We split the y-integral into one over M+ and one over ]R_. Then we can ex- 
change the integrals over M_ and the imaginary axis provided that lie z < fik 
and use dy^vii^k-z) = the same way we integrate over M+ taking 
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z such that < Rez < /Xfc+i. This gives Ay^y^i^k-z) ^ Putting 
these two integrals together we get 



Jr^^ {z - Hk){z - /Xfc+i) ■ ■ ■ {z - Hn)' 



(44) 

which is a constant multiple of e^*^^. This proves (i). For (ii) we proceed 
similarly. Using that 1 < /c, ^ < n we have 



/dx / d. / d^f;^^ (.-...0 



(45) 



When we integrate x over ]R_, we take the z- integral such that Rez < Rety, 
and when we integrate x over M_|_, we choose Re 2; > Re w. Thus, f HSjl reduces 
to 

d^ = 4.. (46) 



27ri /r^,.....^„ (u; - /x^) ■ ■ ■ (w; - /i„) 
Finally, note that 

which shows that the set {^^-k '■ ^ ^ k < n} spans Vn- □ 
Next we verify Assumption (A) from Lemma 3.4 in [6]. Indeed, 

<i^) = 7r^i dw- = CnM<''\^) (48) 

/r^„ - /in 

with c„ = e"*^"/^ 7^ for n = 1, . . . , A^. 

Finally, we can also determine the value of the normalization constant in 
(135|) . since it is given by 1/ det[Mi:i]i<k,i<N with 

Mki = {<Pk*---*<pN* (49) 
Lemma 2.6. VFe /iav e det M = Y[n=i e*'^"'^^, in particular det M > 0. 
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Proof. By Lemma El (i) we may write M^i = {4>k * '^k-i)i^k~^)- Thus, for 
k>i, 

Mk, = f dye^^^y f dze'^'/'-y^ {z - ■ ■ ■ {z - /i,.)- (50) 

Once again, we let run the y-integral over M_ and M+ separately. In the first 
case we take the ^-integral such that Re^; < ^k-, in the second case such that 
Re^; > /ifc. This allows us to exchange the integrals, which gives 

27ri /r^^ z- (Ik 

Since the integrand has no poles for k > i, we have Mke = in this case, 
while for k = i we get Mkk = e"*'^^/^. Thus, M is upper triangular and the 
claim follows. □ 

With the results of Lemma 12. 41 and Lemma [23| Theorem [1] follows directly 
from Lemma 3.4 of [B]. 

We have shown that Theorem [T] holds when we impose /ii < ■ ■ ■ < fi^. 
In particular, the joint density (1201) is given by a A^(A^ + l)/2 correlation 
function: Let m = N{N + l)/2, then 

(|5U]) =m!p(™)({(A5J,n),l < k<n< N}). (52) 

Let M > be any fixed real number. The density fl20p is analytic in each of 
the fij in [— M, M], j = 1, . . . , N. The same holds for the correlation kernel 
(take e.g., = —M — 1). From this it follows that also the r.h.s. of f l52|l 
is analytic in each of the variables /ii, . . . ,fiN- Since this holds for any M, 
by analytic continuation it follows that Theorem [T] holds for any given drift 
vector (/ii, . . . , /i^r) G M^. 



3 Warren's process with drifts 

We have seen in Section [2] that the eigenvalues' density can be written as a 
product of determinants, and, in Lemma [2^61 we calculated the normalization 
constant, so that the probability measure on the eigenvalues reads 



P( n KedAa) =detK'!,(Af)]^<,_,<^n 

^l<fe<n.<Af n=l 



TV 



e 



N 

X 



ndet[0.(Ar\A-)]^^^^^.^„ n dA^ (53) 

n=l l<fc<n<Af 
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Figure 1: The two reflection types in our system. They correspond to the 
boundary condition 

In this section we show that this measure also describes the transition prob- 
abilities for a system of Brownian motions in GT^r, starting from at 
time t = and arriving at x G GTat at time t > 0. Let us consider 
{B^,l < k < n < N} Brownian motions in GTat starting from 0, with 
drift fin, and interacting as follows: 

• the evolution of does not depends on the Brownian motions with 
higher upper index (5™ for m > n + 1, and any i), 

• 5^ is reflected off by 5^-^ and B'^zl 

These reflections are sometimes called oblique reflections [25], since in the 
(x^~^,a;^) plane (resp. {x^Zi,^'^) plane) the reflection direction is not the 
normal one, see Figured] Remark that the projection on {-B", 1 < n < A^} 
differs from the process studied in [21], where the reflections are in the normal 
direction. If we denote by pt the probability density in GT^r of the system 
of reflected Brownian motions, this kind of reflections are described by the 
boundary condition 



whenever = x^"*" or = x^^-^, for 1 < A; < — 1. This is a special 
case of the reflection in a given direction studied in [TH]. For a Brownian 
motion with drift /i reflected at the boundary in the direction v, the boundary 
conditions are the following. Denote by n the normal vector of the boundary, 
let V be normalized such that n ■ v = 1 and let q = v — n. Moreover, set 
D* = n ■ V — q ■ (V — n{n ■ V)) . Then, the boundary condition are written as 



d 



Pt{x) + {fin+l - fJ'n)Pt{x) = 0, 



(54) 



dx^ 



D*'pt = 2(/i ■ n)pt, at the boundary. 



(55) 
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Specializing to our case, we get flM]) . This process, without drifts, was in- 
troduced by Warren in In that paper, he determined the transition 
probabihty for any initial condition and also showed that the process is well- 
defined when starting from 0. We here consider a system of Brownian motions 
with constant (bounded) drifts, which can be expressed as follows, 

Bl{t)=fi,t + bl{t), 

B'^{t)=fLj + b^{t)-L^^-i_^^{t), n = 2,...,N, 

Blit) = /i„t + blit) - L^n-i^^„{t) + L^„„^n-i(t), 2 < A; < n < A^, ^^^^ 

k k k k—1 

Bl{t) = fij + bl{t) + L^r^^^r.-. (t), n = 2, . . . , iV, 

where the b1,l<k<n<N, are independent standard Brownian motions 
and Lx-Y{t) is the local time accumulated at the origin by X{t) — Y{t). The 
question of well-definedness was related with the, a priori possible, presence 
of triple collisions. Bounded drifts do not influence this property as can be 
seen by applying Girsanov's theorem like in the works [T7t [TU] . In this paper 
we do not determine the full transition probability, but rather the transition 
probability for the our specific initial condition. 

Proposition 3.1. Denote by pt : GTn — [0, 1] the transition probability 
density of the N{N + l)/2 Brownian motions with drift, starting from at 
time t = and arriving at x & GTn at time t > 0. Inside GT^v, this density 
satisfies the Fokker-Planck equation 

Moreover, as t \0, the measure pt{x)dx becomes the Dirac measure at 0, 

hmpt{x)dx= Yl ^^k^ (58) 

l<k<n<N 

and the reflections are translated by the boundary condition 

d 

Pt{x) + (/i„+i - fin)pt{x) = (59) 



dx^ 



for 1 < k < N — 1, whenever x^ = x^ or x^ = x^^^ . 

Proof. First observe that by setting ^^'^^(x) = e^^^'^^'^j^lx), we can rewrite 
probability measure on GTat with density 

N N-l n 

pt{x) = det l^^tki^)] Y\e-'^l/' n rre('^"-'^"+iK. (60) 

L J l<fc,£<Af J. J. J. J. 

- - k=l n=l k=l 
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for X = {x^)i<k<n<N £ GTjy. The double product only depends on 
{x^)i<k<n<N-i, while the determinant is a function of {x^)i<k<N- We have 

1 .j,N,t I \ ,j,N,t / \ I 9 ,j,N,t , \ ^N:;i,N,t ( \ 

*7v-fe(a;) = ^ ^7v-fc(a;) + /iw^^7v_fc(a;) - ^^iv-fc(^). (61) 
from which follows that 

(62) 



For /c = l,...,7V — 1, we have 



Pt(x) = (/i„ - /i„+i)pt(x), nX9 Pt(3^) = (/^n. - /in+l) Ptl^^), (63) 



and thus, putting (l62l) and (l63l) together, 
1 ^ " (9^ ^ a 

T\j — 1 Aj — 1 \\i — 1 

\n=l n=l / 

Using that 

TV N~l N-1 N-1 

^/^^~E^"= Xl^(/^n+l-/^n) = ^r^(^n-/in+l)^-2^n/i„(/i„-/i„+i) 
n=l n=l n=l fc=l 

(65) 

the expression between the brackets in fl6^ simplifies to 2 ^ ?7,/i„(/i„ — /x„+i). 
On the other hand, 

N n g N-1 N Q 

^^^^r^^Pt{x) = "^nflnif^n - f^n+l)Pt{x) + flN^Q^Ptix). (66) 
n=l k=l ^ n=l k=l ^ 



Then, (157|) follows from (1641) . ( 165|) and (166|) . The initial condition ( 158|) is 
satified because as t \ 0, we obtain the Dirac measure at x^ = for 
1 < k < N and since we consider pt on GT^r, this immediately implies that 
= for alll</c<n<A^ — 1. Finally, the boundary condition (139]) 
holds trivially by □ 
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n = 5 

n = 4 

re = 3 

re = 2 

re = 1 



^2 



»x|< 



Figure 2: (Left) Initial particles configuration. (Right) A possible particles 
configuration after some time; in this configuration, if particle (1, 3) tries to 
jump, the move is suppressed (blocked by particle (1,2)), while if particle 
(2,2) jumps, then also particles (3,3) and (4,4) move by one unit to the 
right. Particles at level n have a jump rate Vn- 



4 2 + 1 dynamics with different jump rates 

In this section we show that the correlation functions ([3]) that we obtained for 
the GUE matrix diffusion with drifts (and for Warren's process with drifts) 
can be obtained as a limit from an GTjy-extension of TASEP with particle- 
dependent jump rates. This latter process was introduced in [1]. Before we 
come to the convergence result, let us describe the model. 

At a fixed time t, let us denote by x{t) = (x^(t))i<fc<„<Ar G GTjv the 
positions of the A^(A^ + l)/2 particles at time t. We choose initial conditions 
x^(0) = k — n — 1 and let the particles evolve as follows: Each particle has 
an independent exponential clock of rate Vn > 0, i.e., particles on the same 
level have the same jump rates. When the x^-clock rings, the particle jumps 
to the right by one, provided that < x^~^ — 1, otherwise we say that is 
blocked by x^~^. If the x^-particle can jump, we take the largest c > 1 such 
that = x^^l = ■ ■ ■ = x^'^^Zi, and all c particles in this string jump to the 
right by one, see Figure |2] for an example). This ensures that at any time t, 
all the particles are in GTat. More precisely, these dynamics imply that the 
particles stay in a discrete version of GTat, namely 

{{x\ x\...,x^) eZ^ xZ^ X ■■■ xZ^ : x^+^ < xl < (67) 

The joint distribution of the particles has been calculated in Theorem 4.1 
of [5]. It is given by 

N 

const X det [^^/_,{x^)] J] ^")] i<..<n' (68) 

n=l 
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where 



= / d^e*/^z-+^-^(l - v,^,z) ■ ■ ■ (1 - v^z), 

27ri iro (69) 

(f)n{x,y) = ivny~'''lly>^] and y) = 

Actually Theorem 4.1 of |5] is a statement about the marginal of a (possibly 
signed) measure. However, this model is the continuous time limit of a generic 
Markov chain introduced in Section 2 of from which it follows that the 
measure with fully packed initial conditions t/n = x" (0) = -n for 1 < n < iV 
is actually a probability distribution. The form of f l68p follows then from 
the theorem by taking a{t) = t and b{t) = for all t > 0. Also note that 
we put the transition from time t = to time t (which is encoded by 7t,o in 
the theorem) into ^&]v-fc- The correlation functions of this point process are 
determinantal as shown in [S]. What remains to do is the biorthogonalization 
for the generic jump rates. 

Proposition 4.1. Consider a system of particles on GTjv with fully packed 
initial conditions and dynamics described above. Then, at fixed time t, the 
corresponding point process has m-point correlation function given by 

^™((xi, ni), . . . , {xm, rim)) = det[K^{{xi, ra), {xj,nj))]i<ij<rn (70) 
with (xj, rii) G M X {1, . . . , A^} and correlation kernel 



K^{{x,n), {x',n')) = -0(-'"')(x, x') + K'-k{^)<''-ki^% (71) 

where 



k=l 



1 /" "1 71 — Tl 



;n,i 1 / 1 {Z-Vi)---{Z-Vn) 



27ri /r„ ^ (^z - vi) ■ ■ ■ {z - Vk) 



^n,t 1 / 1 w''^'' {w - vi) ■ ■ ■ {w - Ve_i) 



(x) = — dw —— ^ — (74) 

"-^^ ^ 27ri/r^ e*- {w - v^) ■ ■ ■ {w - v^) ^ ' 

Proof. By Proposition 3.1 of [5], we have 

"-'^^ ^ 27ri (1 - i;iti;) ■ ■ ■ (1 - vM ^ ^ 
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for k > 1. A change of variable z = 1/w then yields f l73l) . Next we need to 
verify that {^I'lf, : 1 < k < n} is biorthogonal to : 1 < i < n} (see 

Eq. (3.5) of [5]). We split the sum over Z into two parts, one over a; > and 
one over x < 0. Then, 



x>0 

Jtz ^ ,x+n 



^^J27ri)2/r. fro e'^ z-+-+^ _ . . . (y; _ 

We choose Fq and r^, such that \w\ < \z\ which allows us to put the sum 
inside the integrals. This gives 

= ^Jdw<f ^,tl^il^3!±p±^^. (77) 
(27ri)2 e*- z^ {w - Vk) ■ ■ ■ {w - v^) z-w 

For X < we choose Fq and F^, such that they satisfy \w\ > \z\ which gives 

x<0 

(27ri)2 7r^_^ e*- z^ {w - Vk) ■ ■ ■ {w - v^) z-w 

Thus, 

Finally, we show that {^^i!^ : 1 < I < n} spans the space of functions V^. 
We denote hy ui < ■ • • < the different values of wi, . . . , w„ and the 
multiplicity of Uk, i.e., ai + ■ ■ ■ + = n. Then, we may write 



27ri Jy^ e*"' {w — ni)"i ■ ■ ■ {w — Uy)^^ 



W=Ui \ j-^j^ \ J I 



' 1 d"'"^ /w^+"i-r 

^ (ftj - 1)! dw°»-^ \ p*"" 

i=l ^ ' 

i=l j=l 



^0) 
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For £ = 2, . . . , n, we can represent ^^-e same way, but with exponents 

a£,i < I < i < v. Since (a^,!, . . . , at^u) 7^ ("£,1, • • • , "£,1.) for A; 7^ this 
shows that 

span{$^:^^ ■.l<i<n} = span{x ^ {uiYx^'^^ ■ I <u <v,l < j < aj, 

(81) 

which isVn- □ 

We continue by estabhshing the convergence result under the scahng (j8]). 
Correspondingly, we rescale (and conjugate) the kernel Kt and define the 
rescaled kernel as 

^.",T,esc((e, n), (e', n')) = ^ VtK^- (([rT - ^v^], n), ([rT - ^'v^], n')) 

(82) 

where [■] denotes the integer part, and the drift v is now /i^ = 1 — l^j^jT- 
Of course, T is assumed to be so large that > is satisfied. 

Proposition 4.2. For any fixed L > , the rescaled kernel K^j,^^^^ converges, 
uniformly for ^, ^' G [— L, L], as 

hm i^f.^,_((e,n), (r,n')) = /C"((e,n), (e',n')) (83) 

wt/i ii'^ = Kt- given in (jl]). 

Proof. Let us define the rescaled functions 

^^^.,T,esc(0 = T("-^-+i)/2e-^ + ^Vt), (84) 

<^>:L^.,T,esc(e') = T-("-^)/2e^^ $:Lt(^T + rVT), 

where we also rescale the jump rates as in (|8]). We have to show that 
these functions converge to their analogues from ([5])-(17j). We first verify 
that 0r";e2(^.f) with n < n' . For y > y', the integrand 

of ^'■"''"'^(i/, y') in (172|) has residue at infinity and thus the whole integral 
vanishes, while for y < y', there is no pole at 2 = and therefore 



4.i'-'\y,y') = ± .„(»'-.)+(»'-">-• n-^I„<,1. (85) 
Hence, for its rescaled version, 

n' . . (^_g')VT+(„'_„)_i 

r) = E 1 - ^ 11 ^ iK>eb (se) 
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which, as T — 7- oo, converges to 

1 ( ^\n'-n r p2(^'-0 

y e^»«'-«) TT ^—M,>i'} = / , 

H - /ii 27ri (2; - /i„,+i) ■ ■ ■ (^ - /i„' 

(87) 

Next we show that ^":!fc,T,rcsc(0 ^n-fc(0 uniformly for { G 
We have 

-^^"'^ (88) 
= ^ / dze^^^°(^)+^^^(^)+^^(^)^7(z) 

with fo{z) = 2; — 1 — In 2;, /i(2) = ^ In z, f2{z) = — (n + 1) In z, and 

_ (^(^-i) + ^.)...(Vr(.-i) + ^„)^ 

(v^(z - 1) + fll) ■ ■ ■ (v^(z - 1) + 

A Taylor expansion around the double critical point of /o, i.e., around Zc = 1 
gives 

fo{z) = l{z-ir + Oi{z-lf), 

f,(z)=^(z-l) + 0{{z-lf), (90) 
/2(z) = + 0(^-1). 

Fix r > 1 — -\/T and deform ro,t, to the contour 7 = 71 U 72 with 
7i = 1 -^_/VT + i[-r,r], 72 = {1^1 = r} n {Re z < 1 - ^_/Vt}. (91) 

Let us verify that 7 is a steep descent path for /q. 

On 7i we have Re /o(x + iy) = a; — 1 — | ln(x^ + y^), so that we have 

dRe/o(. + .y) ^^ , ,e[-r.rl, (92) 

dy + 7/^ 

with X = 1 — -\/T. Thus /o is strictly increasing on 1 — + i[— r, 0) and 
strictly decreasing on 1 — + i(0,r]. On 72 we compute 

Re /o(re^''') = r cosy? — 1 — Inr, G (arccos 27r — arccos ^), (93) 

which means that /o is strictly decreasing on 72 fl {Im^; > 0} and strictly 
increasing on 72 fl {Im^; < 0}. Thus 7 is a steepest descent path for /o and 
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the major contribution comes from a line segment 75 = 1 — ^ + i[—6, 6] for 
any 6 G (0, 1). Indeed, the error we make when we integrate along 75 instead 
of 7 is of order 0{e~'^'^) with c ~ 6"^. We therefore consider the integral on 
15 only, 



27ri 



is 



Using |e^ — 1| < |x|e'^', the difference between fl94p and the same integral 
without the error term can be bounded by 



271 



dz 

is 



,ci^VT{z-l)+C2^iz-lf 



xO(^z-l, Vf{z - 1)2, T{z - l)^ T^^-''^/\z - 1)"-'=) (95) 



for some constants Ci and C2 that can be chosen arbitrarily close to 1 as 
5 — >■ 0. By a change of variable Z = \/T{1 — z) one then sees that this error 
is of order 0{T~^^'^). Hence we can consider the integral in flM|) without the 
error term, where the argument of the exponential function simplifies to 



^ / dze-^(^-^)'/2+«-^(^-^)^7(^). 
27ri J^.^ 



(96) 

is 

The error we make if we extend 7^ to 1 — + iM is of order 0{e~'^'^). All 

together the integral from fISS]) agrees, up to an error 0{e~^'^, T"^/^) uniform 
in ,^ G [—L, L], with 

d^e-^(^-i)'/2+«^(^-^)(7(^), (97) 



27ri 



where the poles of g lie on the left of the integration axis. After a change of 
variable Z = —\/T{z — 1) this integral can be identified as '^^-kiO- 
Finally we show that $"l,,r,rcsc(^') ^ K-ki^')- We have 



27ri 



(VT(^; - 1) + ^0 ■ ■ ■ {^{w - 1) + /i,,i) 
(yT(ii; - 1) + //i) ■ ■ ■ ( VT(ti; - 1) + /i,) ' 
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and by a change of variable W = —\/T{w — 1) and a Taylor expansion in the 
exponent we get 

^y-k+l J ^^^^^_^^y2/2+g/H/+0(T-V^) (W - /ii) ■ ■ ■ (W - /Xfc-i) ^ ^gg^ 



2vri Jr^ (w - /ii) ■ ■ ■ (w - /i„) 

which converges uniformly for E,' G [— L, L] to "^'^iT^X'C')- n 

With the above results we can now prove Theorem [3J 
Proof of Theorem\^ Set m = A^(A^ + l)/2 and define ni, . . . , Um by 

ni = 1, = ns = 2, n4 = ns = ne = 3, . . . , nm_Ar+i = • ■ ■ = = A^- 

(100) 

For A C we set Ay = (rT - VTA) n Z. Then, we have 

{xi,...,Xm)&AT 

= ^ d-x det [K^MrT - N Vt, n,), (rT - [x,] v^, n,))] 
= y d"^xdet [^r,r,resc((N,n,),([xj],nj))]^^^^^.^^ 

(101) 

and 



v{A) = J^d^xdet[K;^{{x,,n,), (x,, n,))] (102) 

Since the determinants are continuous functions of the kernels, we have by 
Proposition 14.21 that 

= det[K^{{x,,n,), (x„n,))]^<^^^.<^ (103) 

for all Xi, . . . , Xm £ IR- Thus we have shown that the densities of the proba- 
bility measures in question converge pointwise to each other. Then, ( ITTj) is 
a direct consequence of Scheffe's theorem, see e.g. [3]. □ 
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A Determinantal correlations 

Since we refer several times to this Lemma 3.4 of [6], we report it here. 

Lemma A.l (Lemma 3.4 of [B]). Assume we have a signed measure on 
{x", n=l,...,N,i = l,...,n} given in the form. 



N-l 

— TT det[</.„(a:r,x"+^)]i<,,,<„+idet[^j^_,(xf )]i<,,,<^, (104) 



n=l 



where x^j^i are some "virtual" variables and Zn is a normalization constant. 
If Zn 7^ 0, then the correlation functions are determinantal. 

To write down the kernel we need to introduce some notations. Define 

0(--)(x,y) = l^'^- * ■ ■ ■ * < . (105) 

where (a * 6)(x, y) = J2z& '^(^' ^)K^^ y)> ^''^d, for 1 < n < N , 

^"-,(^):=(0^"'^^*<-,)(2/), J = l,---,iV. (106) 
Set (f)Q{xi,x) = 1. Then the functions 

{(00 * 0(''"^)(a;?,x), . . . , {<Pn-2 * 0("-^'"))(xr?,x),</)„,-i(xr\x)} (107) 

are linearly independent and generate the n-dimensional space Vn- Define a 
set of functions {$"(x), j = 0, . . . , n — 1} spanning Ki defined by the orthog- 
onality relations 

= (108) 

X 

for < i, j < n — 1. 

Under Assumption (A): (j)n{Xn+i,x) = Cn^^Q ~^^\x), for some Cn ^ 0, 
n = 1, . . . , N — 1, the kernel takes the simple form 



"2 



K{m, xi; n2, X2) = -0("-"^)(a:i, Xs) + Kl-ki^i)Kl-ki^2). (109) 



k=l 
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